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We review the relations between the full counting statistics and the field theory of electric circuits. We 
demonstrate that for large conductances the counting statistics is determined by non-trivial saddle-point of 
the field. Coulomb effects in this limit are presented as quantum corrections that can stongly renormalize 
the action at low energies. 



1 Introduction 

The concept of Full Counting Statistics has been introduced in very early days of quantum transport. That 
time, nobody ever thought of such an abstract and complicated problem as the evaluating and measuring 
the higher cumulants of electronic noise. The research has been driven by pure curiosity and has resulted in 
compact and deep Levitov formula. |JJ The significance of this contribution has been underappreciated for 
a number of years. The direct experimental verification seemed to be out of question, while the theorists 
opted to address more primitive and standard problems. 

Nowadays Full Counting Statistics is a reasonably established field attracting attention of many, al- 
though the peak of interest is probably in the past. 12 There are beautiful experiments where the higher 
cumulants have been measured |3| and even single electron transfers have been actually counted. |4| The- 
oretically, the Full Counting Statistics has been evaluated for virtually any important electron transport 
system including even so generic as Anderson impurity model. |5|| 

Still it remains underappreciated that the Full Counting Statistics is eventually not about the marginal 
deviations of electric currents. Let as draw a parallel with general relativity. The general relativity is 
not an art to calculate ridiculously small corrections to Newton's law, although only those can be verified 
experimentally. The general relativity brings us true knowledge about the Universe. Similar to that. Full 
Counting Statistics lies near the heart of quantum theory of electricity and is in fact an indispensible 
element for this. 

The present contribution aims to explicate the link between Full Counting Statistics and quantum field 
theory of electric circuits. The most established example of such theory is the quantum theory of a su- 
perconducting Josephson junction in a dissipative electromagnetic environment ||6l that is readily reduced 
to a single-variable field theory for the superconducting phase across the junction Q. The action of the 
same type governs Coulomb blockade phenomena in non-superconducting systems. ll8ll9l [T0J The physics 
of one-dimensional interacting electrons in the framework of Luttinger model is often reduced to similar 
schemes fTT T?!, where the variable is a drop of a phase over the barrier present in the one-dimensional 
setup and can be associated with the voltage drop at the barrier. 

The structure of the article is as follows. We start (Sections 2,3) by formulating a general quantum 
theory of a simplest electric circuit and see the need and advantage of the FCS in this respect. We show 
that the classical limit of the field theory is not trivial as far as FCS is concerned and obtain the FCS at 
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Fig. 1 a. A simplest electric circuit consisting of two terminals, two connectors A,B and one node, gives rise to a 
quantum field theory for the field in the node. b. Charging energy is represented as a capacitor that provides high- 
frequency cut-off for the field, c. Setup studied in Sections 4,5: Quantum connector (characterized by transmission 
eigenvalues T„) in series with Ohmic one (characterized by the frequency-dependent impedance Z{ui)). 



current bias (Section 4). Quantum effects at big conductances (as compared with the conductance quantum 
Gq = /2ttK) can be incorporated by the renormalization of the action parameters. This is frequently the 
case in field theories fO | We perform the renormalization procedure explicitly in Section 5 for a quantum 
contact in series with an Ohmic one and for a set of quantum contacts connected to a single node (Section 
6). 

The concrete results revisited here were first published in lfT4ll . |fT5l . and |fT6l . 



2 Field theory 

Let us first start with elementary electric circuit theory and reason the quantum extension of it. A circuit is 
made of three archetypal elements: terminals, connectors, and nodes. The voltage is fixed in terminals. A 
connector is characterized by its, I — V characteristics: 

I{t)^I{t-{V{t)}) (1) 

where we have assumed general relation between the voltage and current so that current at the time moment 
t depends on time-dependent voltage at all (previous) time moments, A simplest circuit contains 

two connectors {A and B) in series, so that a single node and two terminals. (Fig. la) Connecting elements 
in this way brings about an extra variable : Voltage Vi (t) in the node. In the elementary circuit theory under 
consideration, this voltage Vi (t) is determined from the current conservation in the node, 

lAit; {Vi{t)}) = lB{t; {V{t) - Vi{t)}). (2) 

assuming the terminal voltages are fixed to and V{t). Once the voltage is determined, one finds the I—V 
characteristics of the whole circuit. Thereby, the full description of the system naturally emerges from the 
two descriptions of the separate connectors. 

Let us try to reason a quantum extension of this theory. First of all, it is convenient to change from volt- 
ages to phases defined as <^{t) = (e /fi)V{t). This allows to treat superconducting and non-superconducting 
systems on equal footing. The phases(voltages) of the terminals can be regarded as time-dependent exter- 
nal parameters while the phase of the node becomes a real quantum variable, that is, an operator 
However, it is hardly an option to formulate the theory in operator formalism since finding the classical 
correspondence becomes a formidable task. Rather, we shall opt for Feynman-Vermon or Keldysh-action 
description of the system where all observables can be presented as path integrals over the time-dependent 
non-operator variable <i>(t). The price to pay for this convenience is quite known: the variable "doubles". 
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The point is that the path integral should performed over two parts of the Keldysh contour that correspond 
to coherent evolution of "kets" and "bra"'s. The variable in principle takes different values {t) , $ j~ (i) 
at these two contours. 

We still want the description of the compound system under consideration to emerge from the two 
descriptions of the separate connectors A and B. What descriptions? Since the variable has "doubled", 
the simple I — V characteristics does not suffice. A connector has to be described by something that 
depends on both variables ^f{t), (t). It is also clear that from all possible quantum variables of the 
connector such description should involve only one: the operator of electric current through the connector, 
/(<), and present the reaction of the current on the variables ^^{t), {t). The proper description of the 
connector is thus provided by the Feynman- Vernon influence functional of two variables = (i) + 

<i>r(i))/2,xw^<i>^w-<frw, 



= ^Tcxp|- I dtmt) + ^xit)\Iit)\pTcxp{- I dt[^^t) + ^xit)\Iit)\). (3) 

Here, the trace is over the electronic degrees of freedom specific for the connector: Thereby they are 
"traced out" and never explicitly enter our field theory. The notation T ( T ) denotes time-ordering of the 
exponentials in ascending (descending) order, these exponents presenting the quantum evolution of density 
matrix p subject to the fields $, x- The easiest way to understand the functional is to expand its log in terms 
of X at X ^ 0. The coefficents of the expansion present the cumulants of the time-ordered current operators 
in the connector that is subject to the external classical phase <i>(t). In particular, lnZ(x = 0) = 0, 
(/(t)) = delnZ/dxit), ((/(ti)/(t2))) = e^d^lTiZ/dxiti)dxit2) and so on. Therefore, the functional 
Z is nothing but the generating function of the current fluctuations. For quantum conductors and slow- 
varying $, X it has been first considered in |[Tl[l7l Now we are ready to build up the quantum description 
of the whole circuit. Since the circuit is nothing but a compound connector, it has to be characterized by 
the similar generating function Za+b that depends on the drop of the fields <i>, x over the circuit. Such 
functional is nothing but a path integral convolution of the functionals of the separate connectors Za and 
Zb, 



^A+s[$,x]= J I?*iPxi^a[*i,Xi]2s[^-*i,X-Xi]- 



(4) 



where the path integration measure I?$il?xi = 11 d.^{t)dx{t)- The overall generating function is the av- 

t 

erage over fluctuating phases $i , Xi at the node of the circuit shared by both conductors. Such convolution 
law is nothing but the presentation of the current conservation in the node. One can see it if one substitutes 
Za+b in the form ^ and carries out the integration over <I>i, xi- 

This is the field theory of the simplest single-node circuit. The extension to a more complicated circuit 
is straightforward. The functional is a product of Z's of all conductors integrated over the extra variables 
$, X defined in each node of the circuit. This functional depends on the phases $ applied in each terminal 
and counting fields x defined in terminals so that it gives statistics of the currents to/from each terminal of 
the circuit. 



3 General Properties and Concrete Connectors 

The difference between classical and quantum effects is not readily manifested in the field theory under 
consideration. Indeed, the current fluctuations in the connectors may be of classical as well as of quantum 
origin but the presentation of these fluctuations is almost the same. For the theory in hand, it is constructive 
to define the difference between classical and quantum as the difference between /ow-frequency and high- 
frequency regimes. This is very much like the theory of frequency-dependent noise: It is known that at 
sufficiently low frequencies kBT/h for equilibrium systems) any noise can be regarded as classical 
irrespective of its origin while quantum mechanics becomes important at higher frequencies. 
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In general the functional dependence of a connector functional Z one the phases may be complicated 
and non-local in time. However, one expects on physical grounds the non-locality to vanish at sufficiently 
low frequencies: The current and its statistics at the time moment t would only depend on the voltage at 
the same moment. Therefore, at sufficiently slow realizations of the fluctuating phases the functional Z 
can be expressed in terms of a single function 5(<i>, x). 



(here we for simplicity specify to non-superconducting systems). 

Let us now concatenate two connectors and try to assess the low-frequency limit of the resulting Za+ b 
taking the path integral in (O. There can be two importantly distinct cases. It may be that the path integral 
is dominated by low-frequency $i, xi so that the result is determined by low-frequency limits of Za-, Zb- 
Moreover, in this case the path integral can be evaluated in the saddle-point approximation with the result 



Here $5 and Xs stand for the (generally complex) values of <i>i and xi at the saddle point where the deriva- 
tives with respect to these phases vanish. Since in field theory the relevance of saddle-point approximation 
is normally associated with classical behavior, we call this classical limit. Generally, high-frequency fluc- 
tuations of the fields also contribute to the path integral. The effect of this contribution is that the actual 

id'] 

Sa+b deviates from the S\^g evaluated with the saddle-point method. Adopting again the common field 
theory terms, we call the deviation the quantum correction. Depending on the connector parameters, the 
correction can be vanishingly small or overwhelming. 

The next statement might look less obvious: The quantum correction present the effect of electron- 
electron interaction in the system and in the context of quantum transport is commonly referred to as 
(dynamical) Coulomb blockade effect. First puzzled question: Where is a capacitor of capacitance C 
providing the charging energy Ec = e^/2C that should accompany any passage about Coulomb blockade? 
Well, it always present (Fig. lb) as a capacitance between the node and ground. The capacitive connector 
contributes to (O with ifTSl 



This term suppresses the high-frequency fluctuations of the fields in the node and is usually needed for 
proper regularization of the theory since it provides the high-frequency cut-off. 

Common wisdom of quantum transport suggests that Coulomb interaction is weak provided the typi- 
cal conductance of the connectors G exceeds by far the conductance quantum and is dominating other- 
wise. One can see this from the estimation of phase fluctuations around the saddle point, (($, x))uj — 
{i{G/GQ)uj + Ecoj'^ /K)'^ . The fluctuation is ^ 1 provided G ^ Gq. We will see that in this regime the 
quantum corrections are small although logarithmically diverge at small frequency. In opposite case our 
system develops a strong Coulomb gap and becomes a sort of SET transistor 

We end this Section with concrete examples of connectors. Ohmic connectors are linear conductors and 
exhibit Gaussian current fluctuations. In terms of Fourier components of the fields. 




(5) 



X) = Sa{^s,Xs) + Sb{^ -^s,X- Xs)- 



(6) 




(7) 



In 2oh: 




{^^{^iujZ-\Lj))x-u. " Xu2ujRc{Z-\uj)) coth{nu;/2kBT)x-^) (8) 



Z{uj) being the frequency-dependent impedance (resistance) of the connector Since the action is of Gaus- 
sian type, the field theory is completely trivial. Any circuits made of Ohmic connectors are reduced to 
Ohmic connectors as well. 
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The Ref. Q has provided the connector action for a tunnel connector between the superconducting 
and/or normal leads, that appeared to be non-Gaussian even for normal case. In the simplest low-energy 
limit such connector is a Josephson junction characterized by Josephson energy Ej and corresponding 
action 

InZj = -i^ J di(cos(2$+(t)) - cos(2$-(t)). (9) 

An arbitrary coherent quantum connector is characterized by the set of transmission coefficients r„. 
The FCS studies UJ have demonstrated that in low-frequency limit the functional corresponds to 

I ± I N 

5A(<i>,x) = |^5(zsgn(<i>))x), 5(e) = 5]ln[l + (e«-l)r„]. (10) 

^'^ n=l 

4 Saddle point and FCS at current bias 

Let us investigate the saddle point of our field theory for a quantum conductor (A) of conductance G in 
series with an Ohmic conductor (B) of low-frequency resistance Z. It is instructive to assume that the 
Ohmic conductor is noiseless, that is, is kept at vanishing temperature. Its contribution to the action thus 
reads (see ^ 

= (11) 

zttZGq 

The subject of interest is how the current fluctuations produced in the quantum conductor are distributed 
in the whole circuit. If the Ohmic resistor is small, ZG <C 1, we have FCS of the quantum conductor at 
voltage bias. Increasing Z to ZG 3> 1, we achieve the current bias for this quantum conductor and access 
the FCS in this limit. 

We apply general Eqs. ^ and (|6]l to our specific circuit that is driven by the voltage source Vq. To 
avoid significant quantum corrections to this field theory (Coulomb blockade effects 0) we assume that 
\Z{uj)\Gq ^ 1 at frequencies Tiui ~ max{eV, ksT) where the quantum corrections start to form. The 
zero-frequency impedance Z can have any value. 

Both the voltage drop V at the quantum conductor and the current / through the conductor fluctuate in 
time for finite Z, with averages / — ^0^(1 + ZG)"^, V = Vo(l + ZG)^^. Voltage bias corresponds to 
ZG ^ 1 and current bias to ZG ^ 1, with lo = Vq/Z the imposed current. There are three characteristic 
time scales: h,/ma.x{eV, kT), e/I, and the i?C-time of the circuit. The low-frequency regime on which 
we concentrate is reached for current and voltage fluctuations that are slow on any of these time scales. 

We seek the cumulant generating function of charge 

Coo \ oo 

j2e'^ip{q)]=Y,m^, (12) 
q=0 ) p=l ^' 

where ((g^)) is the p-th cumulant of the charge transferred during the time interval t. It is directly related 
to the Keldysh action in the saddle point ^ by 

T{i)=TSA+B[eV^|'h.~^i)■ (13) 

To characterize the fluctuations of the voltage across the quantum contact, we will also need the cumulant 
generating function of phase, Q{(,). We use that in the absence of noise in the Ohmic connector V = 
Vf) — ZI. Therefore, Q is related to T{i^ by a change of variables. The relation is 

Q{i) = Y.^mr- = 0oc+^(-^gq$/2), (14) 
p=i 



6 



Yuli V. Nazarov : PCS and FT 



00 being the phase change induced by external vohage in time interval t, <j>o = eV()T/2TTTi. In the 
limit Z — > of voltage bias the saddle point of the Keldysh action is at $i = $, xi = x, and 
from Eqs. (|6]l, (fTZt . and (fl4l i one recovers the results of Ref. [1]: The cumulant generating function 
^oiO — 'T'SAieVo/h, = (j)oS{^) and the corresponding probability distribution 

1 ^ 

n— 1 

We note that the parameter 0o is in fact the number of attempted transmissions per channel. The first few 
cumulants are {q)o = ^oG/Gq, {{q^))o = 0o E„ T^i^ " Tn), {{q^o = 0o E„ T„{1 - T„)(l - 2T„). In 
the single-channel case (N — 1) the distribution ( fTSl l has the binomial form ( l20b . 

After these preparations we are now ready to generalize all of this to finite Z, and in particular to derive 
the dual distribution of phase ( 1211 1 under current bias. Calculating saddle-point values of <i>i, Xi from Eqs. 
© and ( fT3] ) we observe that (z = ZGq) 

m = -[^-^m. <r + zS{a)=^. (16) 

z 

The implicit function cr(^) (which is determined from the saddle point of the action) provides the cumulant 
generating function of charge T for arbitrary series resistance Z. One readily checks that (t>oS{^) 
in the limit z ^ 0, as it should. 

By expanding Eq. ( fTSl l in powers of ^ we obtain a series of relations between the cumulants {{q^Jj 
of charge at Z 7^ and the cumulants {{q^))o at Z — 0. For example, to linear order we find (q) — 
(1 + ZG)~^ {q)o, which nothing but the trivial division of voltage: The mean current / is rescaled by a 
factor 1 + ZG coming from the series resistance. Naively, one may assume that the same rescaling applies 
to the fluctuations. Indeed, to second order one finds ((q^)) = (1 + ZG)^^ {{q^Jjo, in agreement with 
elementary circuit theory. 

However, if we go to higher cumulants we find that other terms appear, which can not be incorporated 
by any rescaling. For example, Eq. ( fTSl l gives for the third cumulant 

^zG 

" (1 + ZGY (1 + ZGf {q)o ■ ^ ^ 

While the first term on the the right-hand-side has the expected scaling form, the second term does not. 
This is generic for p > 3: {{q^)) = (1 + ZG)^'^^^ {{q^)) plus a non-linear (rational) function of lower 
cumulants |fT9l . All terms are of the same order of magnitude in ZG. 

Turning now to the limit ZG ^ 00 of current bias, we see from Eq. (fTSt that T Toe with 

^oo(e) = 90^-905'"" (e A) (18) 

defined in terms of the functional inverse S™^ of S. The parameter go = 0o/^ = Iot/c (which assumed 
to be an integer ^ 1) is the number of charges transferred by the bias current Iq during the time interval 
T. Transforming from charge to phase variables by means of Eq. (O, we find that Q — ^ Qoo with 

GooiO = -qoS"'^-0- (19) 

It is interesting to discuss a single-channel conductor (transmission Ti) separately. In this case the func- 
tional inverse gives the function of a similar form. Eq. ( fT9l l reduces to GociS.) = —Qo ln[l + rj^^(e^'> — 1)], 
corresponding to the Pascal distribution (jJTji. The first three cumulants are {(/>) = qa/Ti, ((0^)) = 
iqo/T?){l Ti), {{4>^)) = (go/Tf)(l - Ti)(2 - T,). 

While the charge Q = qe for voltage bias Vq = h4>Q/eT is known to have the binomial distribution ||T] 

'° 1 T^n _ T, \4>o-q 



^00(9)- ( g (20) 
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we find that the dual distribution of phase $ = 27r(/) for current bias Iq = eqo /t is the Pascal distribution 
lHQl 



- 1 
90-1 



PioW^i: Jr*(i-r)*-«°. (21) 



(Both q and (j> are integers for integer (po and go.) 

In the more general case not depending on the type of the quantum conductor we have found that the 
distributions of charge and phase are related in a remarkably simple fashion for q,(f) ^ oo: 

lnP,(0) =lnP4g)+0(l). (22) 

(The remainder 0(1) equals \n{q/(j)) in the shot-noise limit.) This formula, which valid with logarithmic 
accuracy, is a manifestation of charge-phase duality,|8| and holds for any conductors. 

The binomial distribution (l20l l for voltage bias has the interpretation UJ that electrons hit the barrier 
with frequency eVo/2TTh and are transmitted independently with probability Ti. For current bias the trans- 
mission rate is fixed at lo/e. Deviations due to the probabilistic nature of the transmission process are 
compensated for by an adjustment of the voltage drop over the barrier. If the transmission rate is too low, 
the voltage V{t) rises so that electrons hit the barrier with higher frequency. The number of transmis- 
sion attempts ("trials") in a time t is given by (e/27r?i) Jg V{t)dt = cj). The statistics of the accumulated 
phase (f> is therefore given by the statistics of the number of trials needed for lor/e successful transmission 
events. This stochastic process has the Pascal distribution jTH . 

For the general multi-channel case a simple expression for Pq„ {<j>) can be obtained in the ballistic limit 
(all T„'s close to 1) and in the tunneling limit (all T„'s close to 0). In the ballistic Umit one has Goo{£,) = 
qoi/N + qa{N — g)(e^/^ — 1), corresponding to a Poisson distribution in the discrete variable N(j) — go = 
0, 1,2,.. .. In the tunneling limit Qoo{0 — ~1o — £,/g), corresponding to a chi-square distribution 
Pgo(0) oc (/)*~^e~S'^ in the continuous variable (p > 0. In contrast, the charge distribution P^igiq) is 
Poissonian both in the tunneling limit (in the variable q) and in the ballistic limit (in the variable N(po — q). 

For large go and when the discreteness of these variables can be ignored, we may calculate Pq^ {(/)) 
from Goo{£,) in saddle-point approximation. If we also calculate Pipaiq) from in the same approxi- 

mation (valid for large 00 and q), we find that the two distributions have a remarkably similar form: 

P^,{q) = 7V^„(g)exp[rE(2^0o/T,g/T)], (23) 
PM = Af,o(0)exp[rE(2^0/r,go/T)]. (24) 

The same exponential function 

S(a;, y) = Sa{x, -i^s) - y^s (25) 

appears in both distributions (with the location of the saddle point). The pre-exponential functions 
N^g and N^^ are different, determined by the Gaussian integration around the saddle point. Since these 
two functions vary only algebraically, rather than exponentially, we conclude that Eq. ( l22b holds with the 
remainder 0(1) = ln(g/0) obtained by evaluating lTi[2Tr{d'^T,/dx'^)^/'^{d'^Y,/dy'^)-'^/'^] at a; = 27r0/r, 
y = q/t- 



5 Renormalization by Ohmic connector 

We consider the same circuit and turn to analysis of quantum corrections assuming Z{uj)Gq <C 1 in the 
relevant frequency region. We demonstrate that the main effect of the corrections can be incorporated into 
the renormalization of energy dependence of the transmission eigenvalues of the quantum connector. We 
study this dependence in a non-perturbative limit to obtain an unexpected result: owing to accumulation of 
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quantum corrections, all quantum conductors behave at low energies like either a single or a double tunnel 
junction, which divides them into two broad classes. 

It has been shown that at low energy scales the relevant part of the electron-electron interaction in 
mesoscopic conductors comes from their electromagnetic environment lISTl l9l. The resulting dynamical 
Coulomb blockade has been thoroughly investigated for tunnel junctions flOl. The measure of the inter- 
action strength is the external impedance Z{uj) at the frequency scale il = max(ey, ksT) determined by 
either the voltage V at the conductor or its temperature T. \f z = GqZ{VI) ^ 1 the interaction is weak, 
otherwise Coulomb effects strongly suppress electron transport. 

A tunnel junction is the simplest quantum conductor with all transmission eigenvalues r„ <C 1. Inter- 
action effects for general connectors with T„ ~ 1 are difficult to quantify for arbitrary z. For z <C 1, one 
can employ perturbation theory to first order in z |22|. The contributions Il23ll24l associated the resulting 
interaction correction to the conductance with shot noise properties of the conductor, while the interaction 
correction to noise has been associated with the third cumulant of charge transfer ll25l . To sort this out, 
one shall proceed in the framework of the field theory outlined where all cumulants are incorporated into 
functional dependence of the action on the field x- The recent experiment li26J addresses the correction to 
the conductance at arbitrary transmission. 

A tunnel junction in the presence of an electromagnetic environment exhibits an anomalous power-law 
I-V characteristic, I{V) ~ The same power law behavior is typical for tunnel contacts between 

one-dimensional interacting electron systems, the so-called Luttinger liquids lill|. It has also been found 
for contacts with arbitrary transmission between single-channel conductors in the limit of weak interactions 
ll27l . In this case, the interactions have been found to renormalize the transmission. 

In our model of a quantum connector, its transmission probabilities r„ are energy independent in the 
absence of interactions. We first analyze the quantum correction to first order in z. We identify an elas- 
tic and an inelastic contribution. The elastic contribution comes with a logarithmic factor that diverges 
at low energies suggesting that even weak interactions can suppress electron transport at sufficiently low 
energies. To quantify this we sum up quantum corrections to the action in all orders in z by a renormal- 
ization group analysis. We show that the result is best understood as a renormalization of the transmission 
eigenvalues similar to that proposed in lIZTl . The renormalization brings about an energy dependence of 
the transmission eigenvalues according to the flow equation 



To calculate transport properties in the presence of interactions, one evaluates Tn{E) at the energy £^ ~ il. 

With relation (|26] | we explore the effect of quantum corrections on the distributions of transmission 
probabilities for various types of mesoscopic conductors. In general, their conductance G and their noise 
properties display a complicated behavior at z| ln£'| ~ 1 that depends on details of the conductor How- 
ever, in the limit of very low energies z\ ln£'| ^ 1 we find only two possible scenarios. The first one 
is that the conductor behaves like a single tunnel junction with G{V) ~ V"^^ . In the other scenario, the 
transmission distribution approaches that of a symmetric double tunnel junction. The conductance scales 
then as G{V) ~ V^- . Any given conductor follows one of the two scenarios. This divides all mesoscopic 
conductors into two broad classes. 

We still analyze a simple circuit that consists of a mesoscopic conductor in series with an external 
resistor Z{lo) biased with a slow-varying voltage source Vo(t) (Fig. 1) but now concentrate on quantum 
corrections. 

As we have already done, we present the generating function Z([x, $]) of the low-frequency current 
fluctuations in the circuit is as a path integral over the fields $i (i) , Xi (^)(Eq. IHi. It is convenient for us to 
change the order of the connectors so that 



dT.n{E) 



2zTr,{E)[l-T,,{E)]. 



(26) 



(ME 




(27) 
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where d^{t) /dt = eV{t) and Zohm is given by (O. Let us assume that $, x are slow fields. In order not 
to repeat the considerations of the previous Section, we will simply set Z{{)) to 0. In this case, the saddle 
point is trivial: $i, Xi — 0. In physical terms, all the voltage drops on the quantum contact. 

We start the renormalization by concentrating on the "fast" part of the fields <i>i,xi and expanding 
the action till quadratic terms in these fields. Doing so, we neglect the time-dependence of slow fields in 
comparison with that of fast fields, so the corresponding part of the action reads 



(28) 



where a, (3 = ij^ij = ^li^o) ± Xi(^)/2, and Mohm presents (O. We require that Z{uj)G <C 1 at 
any frequency. Under these conditions, the fluctuations of the fast fields are determined by the Ohmic 
term while the part of the action that comes from the fluctuation and does depend on the slow fields $, x is 
determined by the quantum conductor Indeed, taking the Gaussian integral (l28T l we obtain the contribution 
to the action 



<55e($(t),xW) 



dujS^S,mt),x{t)) 



2-K 



a/3 



27r ?,A,°'f)fh^ 



[Nr^yf^ {uj) (29) 



This is the renormalization sought. The correction to the conductance of the quantum conductor it gives is 
of the order of Gz. 

To proceed, we need the action of the quantum conductor at fast fields, not just at slow ones as given 
by Eq. [lOl It is expressed in terms of Keldysh Green functions Gr,l (the "check" denotes 2x2 matrices 
in Keldysh space) of electrons in the two reservoirs adjacent to the conductor ifTTl . It takes the form of a 
trace over frequency and Keldysh indices. 



1 + |({Gl,Gh}-2) 



(30) 



and depends on the set of transmission eigenvalues T„ that characterizes the conductor. The fields 0* {t) 
enter the expression as a gauge transform of G in one of the reservoirs. 



Gr^G"'' and GL(t,i' 
e#^(t) 







i0c (*) 



G'°'(t-i') 



e^"^-'^^"^ 



(31) 



(jx^ — $±x/2 — cf)^™ being the drop of the phase over the quantum conductor, G'"^^ being the equilibrium 
Keldysh Green function 



G^^'^ie) 



1 - 2/(e) 2/(e) 
2[l-/(e)] 2/(e)-l 



(32) 



at a given equilibrium electron distribution function /(e). 

To zeroth order in z the fields cjx^ (t) do not fluctuate and are fixed to eVt ± x/2. Substituting this into 
Eq. (|30] | we recover the slow-field action (fTOl l 



(33) 
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(ffj = f and fi (e) = /(e ~ eV)). To assess the renormalization correction, we expand the non-linear 5c to 
second order in the fluctuating fields (p^ and use ( |29] l. The expression for the correction can be presented 

as 

S^'\V, x) = dio [[2N{oj) + \]S^^ + N{oj)S^;>{i,) + [N{l,) + (34) 

The three terms in square brackets correspond to elastic electron transfer, inelastic transfer with absorption 
of energy hu) from the environment, and inelastic electron transfer with emission of this energy respec- 
tively. It is crucial to note that inelastic processes can only occur at frequencies lo < Vl and that their 
contribution to the integral is thus restricted to this frequency range. In contrast, elastic contributions come 
primarily from frequencies exceeding the scale O. If z = const(<j;) for cj < A, the elastic correction di- 
verges logarithmically, its magnitude being ~ z In A/fi. This suggests that i. the elastic correction is more 
important than the inelastic one and ii. a small value of z can be compensated for by a large logarithm, 
indicating the breakdown of perturbation theory. The upper cut-off energy A is set either by the inverse 
i?C-time of the environment circuit or the Thouless energy of the electrons in the mesoscopic conductor 
The concrete expression for s[^^ reads 

5«(^,X)=*E / ^A.i?+{r„(/L-/+) + 2T„(e^'^-l)/i(l-/+) 

+2T2(cosx - - iWi - Jr) + TnDn + (1 - - At)} 

+ {R^L,x^-x}, (35) 

where we have introduced the functions 

Dn = {l + Tn [/l(1 - IrW - 1) + - /L)(e-^'^ - 1)] (36) 
and the notation 



/+(£) = /(£ + a;), D+{e) = D,,{e + Lo). (37) 

We do not analyze further and instead turn to the analysis of the elastic correction. It is important that 
the explicit form of this correction can be presented as 

^ci" = ^^^^r ^^^^ " "^^"^^ " ^^^^ 

n ^ 

This suggests that the main effect of renormalization is to change the transmission coefficients Tn- It also 
suggests that we can go beyond perturbation theory by a renormalization group analysis that involves the T„ 
only. In such an analysis one concentrates at each renormalization step on the "fast" components of (fr^ with 
frequencies in a narrow interval 5u) around the running cut-off frequency E. Integrating out these fields one 
obtains a new action for the slow fields. Subsequently one reduces E by Slo and repeats the procedure until 
the running cut-off approaches 57. We find that at each step of renormalization the action indeed retains 
the form given by Eq. ( |30] l and only the r„ change, provided z <C min{l, Gq/G}. The resulting energy 
dependence of the Tn obeys Eq. ( |26] l. The approximations that we make in this renormalization procedure 
amount to a summation of the leading logarithms in every order of the perturbation series. 

In the rest of the Section we analyze the consequences of Eq. ( |26] | for various mesoscopic conductors. 
Equation (|26] l can be explicitly integrated to obtain 
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in terms of the "high energy" (non-interacting) transmission eigenvalues T . A mesoscopic conductor 
containing many transport channels is most conveniently characterized by the distribution pa{T) of its 
transmission eigenvalues |28|. It follows from Eq. ( [39l ) that the effective transmission distribution at the 
energy scale E reads 



We now analyze its low energy limit ^ ^ 0. Any given transmission eigenvalue will approach zero in 
this limit. Seemingly this implies that for any conductor the transmission distribution approaches that of 
a tunnel junction, so that all r„ <C 1. The overall conductance would be proportional to ^ in accordance 
withRef. 1221. 

Indeed, this is one of the possible scenarios. A remarkable exception is the case that the non-interacting 
PA has an inverse square -root singularity at T ^ 1. Many mesoscopic conductors display this feature, 
most importantly diffusive ones ll28ll . In this case, the low-energy transmission distribution approaches a 
limiting function 



The conductance scales like is known to be the transmission distribution of a double tunnel 

junction: two identical tunnel junctions in series (29l. Indeed, one checks that for a double tunnel junction 
the form of the transmission distribution is unaffected by interactions. This sets an alternative low-energy 
scenario. We are not aware of transmission distributions that would give rise to other scenarios. 

We believe that this is an important general result in the theory of quantum transport and suggest now a 
qualitative explanation. The statement is that the conductance of a phase-coherent conductor at low voltage 
and temperature f2 ^ A asymptotically obeys a power law with an exponent that generically takes two 
values. 



For tunneling electrons the exponent is 2z. An electron traverses the conductor in a single leap. The second 
possible exponent z has been discussed in the literature as well, in connection with resonant tunneling 
through a double tunnel barrier in the presence of interactions [ TT^. This resonant tunneling takes place via 
intermediate discrete states contained between the two tunnel barriers. The halved exponent a ~ z occurs 
in the regime of the so-called successive electron tunneling. In this case, the electron first jumps over one 
of the barriers ending up in a discrete state. Only in a second jump over the second barrier the charge 
transfer is completed. Since it takes two jumps to transfer a charge, the electron feels only half the counter 
voltage due to interactions with electrons in the environmental impedance Z at each hop. Consequently, 
the exponent at each jump takes half the value for direct tunneling. Our results strongly suggest that 
this transport mechanism is not restricted to resonant tunneling systems, or, in other words, that resonant 
tunneling can occur in systems of a more generic nature than generally believed. As far as transport is 
concerned, a mesoscopic conductor is characterized by its scattering matrix regardless of the details of its 
inner structure. In this approach it is not even obvious that the conductor can accommodate discrete states. 
Nevertheless, the transmission distribution of this scattering matrix does depend on the internal structure of 
the conductor The inverse square root singularity of this distribution at T ^ 1 for a double tunnel barrier is 
due to the formation of Fabry-Perot resonances between the two barriers. Probably similar resonances are 
at the origin of the same singularity for more complicated mesoscopic conductors with multiple scattering. 
They are then the intermediate discrete states that give rise to the modified scaling of the conductance in 




(40) 




(41) 




(42) 
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presence of interactions. One may speculate that in diffusive conductors these resonances are the so-called 
"nearly localized states" found in | 30l. 

From equation (3% one concludes that the resonant tunneling scaling holds only if G{E) ^ Gq so 
that many transport channels contribute to the conductance. At sufficiently small energies, G{E) becomes 
of the order of Gq. All transmission eigenvalues are then small and the conductance crosses over to the 
tunneling scaling. 

6 Renormalization by quantum connectors 

One could wonder about the generality of the results obtained in the previous Section and expressed by 
Eq. |26]lndeed, it has been proven under rather restrictive assumptions of an Ohmic connector of negligible 
resistance. Here, we consider a more general model of several quantum connectors coming together in a 
single node (Fig. |2]i. Each connector labeled by k is characterized by the set of the transmission eigen- 
values Tn ■ In traditional Coulomb blockade situation (G <C Gq) this setup is called Coulomb island or 
SET transistor |9| and is seen very different from the junction-in-the environment setup considered in the 
previous Section. However, we show that in the limit of large conductance G ^ Gq the Coulomb island 
is governed by very similar renormahzation equations: 



It looks like each quantum channel sees all others as an "environment" characterized by the effective island 
conductance g = J2n k Tn^- If we use z — g^^, Eqs. I43]andl26lare identical. The difference is that g by 
itself is subject to renormalization. 

This results in very different low-energy behavior In contrast to the considerations of the previous 
Section, the renormalization of all transmission eigenvalues may break down at finite energy — effective 
Coulomb gap — Ec oc goEce~°'^° , a being a numerical factor depending on the details of the initial trans- 
mission distribution, go is the island conductance at high energies > Ecgo- Remarkably, Ec coincides 
with the effective charging energy evaluated with instanton technique. ||3T1 However, the renormalization 
stops at the effective Thouless energy ETh ^ G{E)6 /Gq, 6 being mean level spacing in the island. This 
gives rise to two distinct scenarios at low energy. If go > a^^ \n{Ec/S), Coulomb blockade does not 
occur with zero-bias conductance being saturated at the value G(£^Th) ^ Gq. Alternatively, G(0) « 
and Ec defines the Coulomb gap. 

Let us give the details of the model in use. The Coulomb island is characterized by two parameters: 
charging energy Ec and the mean level spacing S, Ec ^ S. Il32l The island is connected to M > 2 ex- 
ternal leads by means of M arbitrary quantum connectors (Fig. |2]l characterized by the set of transmission 
eigenvalues Tn^ We assume that the island is strongly coupled to the leads, go = J2n m ^- ^™ 

goal is to evaluate the functional Z for the whole circuit that now depends on voltages and counting fields 
in each terminal, 2{[Vi, Xi]))- To evaluate this for the Coulomb island, we have extended the semiclassical 
approach for the FCS of the non-interacting electrons f33 |. The node houses a dynamical phase variable 
(j){t) 1 18 1 , its time derivative, (f){t) / e, presents the fluctuating electrostatic potential of the island. Accord- 
ing to the rules of our field theory, the functional is represented in the form of a real-time path integral over 
the fields (j)^ {t) residing at two branches of the Keldysh contour 




2Tr'(i-rro 



(43) 



dlnE 




— oo 




(44) 



fc 
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Fig. 2 Coulomb island setup: A single node connected to M = 3 terminals by quantum connectors. The field theory 
is for the fields 0* (t) . 



Here $ = [ J'^'' , ^, . ) is the matrix in Keldysh space, 2x2 matrix G(ii, t2) presents the electron 





Green function in the island that implicitly depends on (f>^{t). The trace operation includes the summation 
over Keldysh indices and the integration in time. The contribution of each connector Smn has a form (|30] | 



Sil = -i 5] Tr In [l + i7;W ({G, Gf} - 2) 



(45) 



{G, G^} denoting the anticommutator of the Green functions with respect to both Keldysh and time in- 
dices. The Green functions in the leads Gk{x) obtained by x-dependent gauge transformation |34| of 
the equilibrium Green functions in the reservoir k, g|°', G^(e) — e:xjp{ixk f3/2)G^^'' (e) exp(— jxfcf3/2), 

where gI*^ are given by g{,"' = ^ ^(-^ ) 2/^^'^l "^^^"^ /fc (^) P''^^^"'^^ ^^^^ ^1^'-^°" '^i^^i^'^tion 
function in the fc-th reservoir. The expression ( |45] ) is valid under assumption of instantaneous electron 
transfer via a connector, thus corresponding to energy-independent Tn ■ 

In order to find G(<i, at given {t), we minimize the action with respect to all G{ti, subject to 
the constrain G o G = 5{ti ~ 12). This yields the saddle point equation for G[ti,t2)'- 

tt^ + Tl'^ ({Gl,G}-2] 



where [..,..] denotes the commutator in the Keldysh-time space. This relation expresses G(ti,t2) = 
G(ti, 12] [4>^{t)]) via the reservoir Green functions This circuit theory relation is similar to obtained 
in ||33J . It disregards the mesoscopic fluctuations, since those lead to corrections of the order of l/go 
at all energies, whereas the interaction corrections are of the order of ^ l/gg In(-E) tending to diverge at 
small energies. If (t>±{t) = 0, Eq. (l46T l separates in energy representation and coincides with that of Ref. 
ISl. 

This sets the model. We start the analysis of the model with perturbation theory in (j)± around the 
semiclassical saddle point G{ti,t2) = Go, 4'±{t) = 0. The phase fluctuations are small, (50^ ~ l/so^ 
so we keep only quadratic terms to the action (l44l i. The resulting Gaussian path integral over (f)± can be 
readily done. This procedure is equivalent to the summation of all one-loop diagrams of the conventional 
perturbation theory, i.e. to the "random-phase approximation" (RPA). 
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We restrict ourselves to the most interesting low voltage/temperature limit, max{eV, kT} ^ goEc- In 
this limit, we evaluate the interaction correction to the CGF with the logarithmic accuracy. It reads 

go \max{eV,kT} 
de ^ 2tJ/-1(1 - tI!'^){{GI, Go} - 2) 



AS-^ = — In ( r^,^ ,,^1 ) X (47) 



2- ft 4 + i;['=l({G^,Go}-2) 

provided max{ey, fcT} > i?Th, where i?Th — goS is the Thouless energy of the island. In the opposite 
case, ma.x{eV, kT} < Etii, the voltage/temperature should be replaced with i?Th- Note, that the cor- 
rection ( l47b is contributed by only virtual inelastic processes that change the probabilities of real elastic 
scatterings. 

For simplicity, we consider the shot-noise limit eV 3> kT only. Then the magnitude of the correction 
shall be compared with the zero-order CGF 5"^°' ~ toeVgo. This implies that the perturbative RPA 
result ( |47] l is applicable only if ^q"^ In {g^Ec / eV) <C 1. At lower voltages logarithmically diverges. 
This indicates that we should proceed with a renormalization group (RG) analysis. 

We perform the RG analysis of the action ( |44] | along the lines of the previous section decomposing 
(fy^ [t) onto the fast and slow parts . On each step of RG procedure we eliminate the fast degrees of 
freedom in the energy range E — SE < u; < E to obtain new action Se-se[4>s], E being the current 
ultraviolet cutoff. Our key result is that the change in the action at each step of RG procedure can be 
presented as a change of transmission eigenvalues TA ' . Therefore, the RG equations can be written directly 
for transmission eigenvalues and take a simple form (03]). The equations are to be solved with initial 
conditions at the upper cutoff energy E = goEc, those are given by "bare" transmission eigenvalues 
Tn'^ {E — goEc) — Tn^ . The RG equations resemble those for the transmission coefficient for a scatterer 



in the weakly interacting one-dimensional electron gas 112711 and for a single multi-channel scatterer in the 
electromagnetic environment lfT5l . The effective impedance Z is just replaced by inverse conductance of 
the island to all reservoirs, G{E) — Gq J2n k '^'^^ (^)- The important difference is that this conductance 
is itself subject to renormalization. The difference becomes most evident in the case when all contacts are 
tunnel junctions, tIi'^ ^ 1. In this case, one can sum up over fc, n in Eqs. ( |43] l to obtain the RG for the 
conductance only : dG/d\nE = 2Gq. This renormalization law f35l was recently applied to conductance 
of granular metals. The Eqs. ( |43] | could be also derived in the framework of functional RG approach to 
(7-model of disordered metal. Il36l . 

We solve the RG Eqs. ( |43] ) in general case to obtain 

rW(i?) = rWy/(l-7l'=l(l-y)), (48) 
\nigoEc/E) = -i ^ ln(l - (1 - y)) (49) 

n, k 

The first equation gives the renormalized transmission eigenvalues at a given value E of the upper cutoff 
in terms of variable y{E), < y < 1. The second equation implicitly expresses y{E). 

We note that the energy dependence of transmission coefficients induced by interaction is very weak 
provided G{E) ^ Gq: If energy is changed by a factor of two, the conductance is changed by ^ Gq. 

To use the equations for evaluation of FCS at given voltages y^'^l of the leads, one takes tI^^ (E) at upper 
cutoff E = maxj: ( V^t*^! ), and further disregards their energy dependence. Then one can follow the lines of 
Ref. lf33l : It is convenient to introduce the function (x) = — J2n ln[l + ^Tn^ {x — 1)] to incorporate 
all required information about transmission eigenvalues. The renormalization of S'^'^'l in terms of y is 
especially simple: S^''^{x,y) = S^''\{x + l)y - 1) - S^''\2y - 1). From this one readily finds the 
conductance of each scatterer, G^^\y) = 2GQdS^''^ /dx{l,y), as well as the renormalized transmission 
distribution TV''"' (T, y) = (2/7r)Im{aS'W /dx (1 - 2/T - iO, y)}. 
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Energy E/E^ \ Energy E/Ec 



(a) (b) 

Fig. 3 The total conductance of the Coulomb island versus the energy: two scenarios. We assume \n{Ec/S) — 10.0. 
Arrows show the energy scale ~ S. Pane (a): tunnel connectors, go changes from 42 (upper curve) to 14 (lowermost 
curve) with the step 4. Pane (b): diffusive connectors, go changes from 18 to 6 with the step 2. The conductance either 
hits manifesting the Coulomb gap or saturates at finite value. 

\k] 

The RG equations (1) have a fixed point at TA = 0, y = that occur at finite energy 

E^Ec^ goEc Y[il ^ T W ) (50) 

This indicates the breakdown of RG and formation of Coulomb blockade with the exponentially small 
gap Ec- The same energy scale was obtained from equilibrium instanton calculation of Ref. (31). For a 
field theory, one generally expects different physics and different energy scales for instantons and pertur- 
bative RG. The fact that these scales are the same shows a hidden symmetry of the model which is yet to 
understand. 

Alternative low-energy behavior is realized if the current cut-off reaches _ETh — G{E)S /Gq. (Fig. 
3) The log renormalization of the transmission eigenvalues stops at this point and their values saturate. 
We thus predict a sharp crossover between the two alternative scenarios, that occur at value of go = 9c 
coiTesponding to Ec — 6. This value equals gc — a^^ hi{Ec/5), where a = ^gQ^ J2n k ^^(^ ~ ^n'^')' 
and depends on transmission distribution of all connectors. If all connectors are tunnel junctions, ar = 2. 
For diffusive connectors, an — tt^ /8 and the energy dependence of the total conductivity is given by 
9d{V) ^ goV^ctgV?, e = 250-1 HgoEc/eV). (Fig. 3) 
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